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ABSTRACT 

We present a study of the luminosity density distribution of the Galactic bar using 
number counts of re d clump giants (RC Gs) from the OGLE-III survey. The data were 
recently released by iNataf et al. I (120121) for 9019 fields towards the bulge and have 



2.94 x 10 6 RC stars over a viewing area of 90.25 dcg 2 . The data include the number 
counts, mean distance modulus (/x), dispersion in \x and full error matrix. We fit 
the data with several tri-axial parametric models, after excluding the bad sightlines, 
2.65 x 10 6 RC stars arc remaining with the sky coverage of 81.47 deg 2 . We use the 
Markov Chain Monte Carlo (MCMC) method to explore the parameter space and find 
that the best-fit model is close to exponential (E2 model), with the distance to the 
GC is 8.17 kpc, the ratio of semi-major and semi-minor bar axis scale lengths in the 
Galactic plane xq, yo, and vertical bar scale length zq, is xq : yo ■ zq w 1.00 : 0.41 : 0.38 
(close to being prolate). The bar major axis is ~ 0.68 kpc and has an angle of between 
29° and 32°, slightly larger than the value obtained from a similar study based on 
OGLE-II data. The number of estimated RC stars is 2.51 x 10 6 , which is systematically 
lower than the observed value. We subtract the smooth parametric model from the 
observed counts and find that the residuals are consistent with the presence of an 
X-shaped structure in the Galactic centre, the excess to the estimated mass content is 
~ 5.3%. The results can be used as a key ingredient to construct new density models 
of the Milky Way and will have implications on the predictions of the optical depth to 
gravitational microlensing and the patterns of hydrodynamical gas flow in the Milky 
Way. 

Key words: Galaxy - structure: Galaxy - stellar content - Galaxy: bulge - Galaxy: 
centre 



1 INTRODUCTION 



Most spiral galaxies in the local Universe are barred. Galactic bars are often straight and their presence appears to be 
intimately related to star formation. Rings are often seen in these galaxies and they are likely related to resonances in a 
rotating barred potenti al. Thus barred galaxies offer exciting opportun ities to understand gas and dynamical processes in 
spiral galaxies (see e.g. Sellwood fc Wilkinson 19931 : Athanassoula 20121 for reviews). It is now generally accepted that the 
inner region of our own Galactic also hosts a tri-axial, b ar-like structure. O bservational evidence for a bar has arisen from 

Binnev et al. 1991 ). s urface brightness (e. g. Blitz fc Spergell 



several sources, s uch as the studies of gas kinematics (e 
star counts (e.g. iNakada et al.lll99ll ; IStanek et al.lll994l ) and microlensing (e.g. lUdalski et all 11994 : see 

review' 



Gerhar 



d 



1991 



2002 for a 



Various tracers have been used to constrain stellar density models of the inner G alaxy. iDwek et al. I (|l995h used the COBE- 
DIRBE multi-w avelength observatio ns of the Galactic centre (jWeiland et al.ll 19941 ) to constrain several multi-parameter ana- 
lytic bar models. Stanek et al. ( 19971 ) first used red clump giants (RCGs) to constrain bar models (see als o Babusi aux fc Gilmore 
20051 and iNishivama et al.112005 who traced the bulge RCG population in the infrared). iRattenburv et al] (|2007l ) applied the 
method of Stanek et al. ( 1997t ) to OGLE-II data and updated the bar parameters using the RCG population. 
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Unfortunately these studies leave many bar parameters uncertain. This may reflect the intrinsic co mplexities of the 



Galac tic bar. For example, there has been a fair amount of discussion about the existence of a nuclear bar. iNishivama et al 



( 2005h bund a clear change of slope in the red clump giants around \l\ = 4°. which was confirmed recently by the VVV 



survey (jGonzalez et al. II2011T ). This was interpreted as the existence of a nuclear ba r. However, numerically simulated bar s 



appear to be able to reproduce the changes in the slope without any nuclear structure (jGerhard fc Mart inez- Valpucsta 1 l2012l ) 



With more and mor e dynamical dat a emerging from seve ral kinematic surveys (e.g., BRAVA, 



Howard et al 



2008, ARGOS, 



Freeman et al. 20131 . and APOGEE, Nidever et al. 2012 ). it appears increasingly important to construct a more accurate 



photo metric model, which will serve as a key ingredient in self-consistent dynamical modelling (|Wang et al.ll2012l ; lLong et al 
20131) . 



The outline of this paper is as follows. In iJ5]we summarize the observational data used, as released bv lNataf et al.l (|2012f ); 
in §3] we describe the methodology used in this analysis; in ij4] we apply our methods to the OGLE-HI data, and present the 
main results of this work. In ij5]we discuss our results and compare them with earlier work. 



2 DATA SAMPLE 



In this paper, we explore the OGLE-III observation data recently released by iNataf et alj (120121 1 who carefully studied the 
extinction and the red clump giant populatio n using both OGLE and 2MASS data. T he RC is a prominent, well-populated, 
and localized feature of Galactic bulge CMDs (jTerndrup et al. 1998; Stanek et al ]| 19941 ). RCGs are very good standard candles 
with Mi — —0.12 and an in trinsic dispersion of 0.09 mag in luminosity at the estimated age and metallicity distribution of 
the bulge stellar population! Nataf et al. 20121). 



The final released data of lNataf et alj (|2012D contain 9019 sightlines, which are almost entirely within the range —10° < 
I < 10° and 2° < |6| < 7° and have a sky coverage of 90.25 deg 2 . The released catalogue contains the total number counts, 
mean distance modulus, geometrical variance, i.e. square of the dispersion in distance modulus, and full error matrix for each 
sightline. 



These new data address two shortcomings in the study bv lRattenburv et al.l (|2007l ). First, the total area in lRattenburv et al 

1 20071 ) covered only about 11 square degrees, spread over a longitude range of —6.8° ^ I ^ 10.6° degrees. Furthermore, the 



field la titudes are clustered around (—4°) with only three fields at positive latitudes. Second, more seriously. iRattenburv et al 



1 20071 ) found that in order to m atch the luminosity of local and bulge red clump giants, an offset of 0.3 magnitude had to be 
artificially applied. The study by Nataf et al. ( 20121 ) addressed both issues. First, they resolve the "offset" difficulty by a much 
more careful consideration of extinction and a new calibration in the intrinsic luminosity of the red clump giants. Second, 
the number of fields is significantly m ore (9019), covering larger than 90 square degrees, nearly one order of magnitude larger 
that used in IRattenburv et alj (|2007l ). As a result, the differences between different models become much more drastic and 
the constraints on the model parameters are tighter. The data also allow us to explore the residuals between data and model 
much more convincingly, which we will exploit in some detail (see E0J) - 



3 METHODOLOGY 

In this section, we first describe the parametric tri-axial models we use to model the luminosity density model of the Galactic 
bar, then we outline the Markov Chain Monte Carlo method we use to find the best-fit parameters for each model. The results 
of our procedure will be presented in the next section. 



3.1 Parametric models 

We first re-examine a few models of bar density distribution proposed by Dwek et al. 1 19951 ): 
nd = N exp(-r 2 /2), 
riE 2 = N exp(-r), 
riE 3 = N K (r B ), 

where Ko is the modified Bessel function of the second kind and 

2 -1 V2 



1/4 



(-) + " 

V xo J \ yo I \zo 

(.-)"■ 



(1) 

(2) 
(3) 

(4) 
(5) 
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These models are representative of a larger suite of parametric models studied by Rattenburv et al. ( 20071 ) . 

We set the origin of the coordinate system at the Galactic Centre (GC), with the x-y plane aligned with the mid-plane of 
the Galaxy and the z-direction parallel to the direction of the Galactic poles. The x-direction is defined to be co-linear with 
the semi-major axis of the bar. The functions are rotated by an angle a around the z-axis. iVo is the number counts per kpc 3 , 
and the distance to the GC is another free parameter _Ro- Thus we have six parameters to model the bar structure. 

The bar can also be rotated by a tilt angle /3 around the y axis, corresponding to the Sun's position away from the 
mid-plane of the Galaxy. Furthermore, there is tentative evidence that the centroid of the bar may be offset from the centre 
of mass of the Galaxy (as evidenced, for example, from the possible gamma-ray line emission offset from the centre, see 
Su fc Fmkbeinerll2012r i. We thus include two additional parameters in the modelling process, SI, 8b, which allow for this 



possible centroid offset. We also fit the data using models in which all nine parameters are allowed to vary and compare the 
results with those found using the simpler six-parameter models. 

With these models, we can calculate analy tically the number counts, mean distance modulus and geometrical variance of 
each sightline as derived bv lNataf et alj ( 20121 1. 

The number counts in each sightline Ni is: 



n(s)s ds, 



(6) 



where the constant a depends on the solid angle subtended around each sightline, and s m i n and s max are the minimum and 
maximum distances along the line of sight. We take s m in = 3 kpc and s max = 13 kpc. We perform the integration over the 
range of ~ Ro ±5kpc as we do not expect the tri-axial bulge density structure to exceed these limits. The results are insensitive 
to the choice of integration limits since the bar density falls away rapidly and the contribution to the integral over the outer 
part of the bulge is negligible. Taking ~ Ro ± 3kpc instead, the total number count changes only by 0.2% on average, and 
thus has the integration limit has little effect on our results. 

The mean distance modulus of each sightline /i; can be calculated as: 



n(s)s fids , 



n(s)s 2 ds, 



where the distance modulus is n = 51og(s/kpc) + 10. 
The geometrical variance cr Ml j can be written as, 



n(s)s fj, ds , 



n{s)s ds 



Mi 



(7) 



We can evaluate the \ 2 goodness-of-fit statistic using the number counts, mean distance modulus, geometrical variance 
and their covariance matrix. The \ 2 of the bar model can be written as, 

7-obs \ 



X 



= ^(xr d - x° bs ) T sr i (xr d - x° 



(9) 



Xf od comprises the number counts (Ni), mean distance modulus (in) and geometrical variance (crjj;) for the i-th sightline 
for the model under consideration, while X° bs comprises the observed values of these variables for the s ame sightline. S; is 
the covariance matrix of uncertainties. The data comprising X™ od are available from Nataf et al. |2012r i. The data include 
9019 lines of sight in total, but not all are used, we exclude the "bad" sightlines which deviate far from the model (see below). 



3.2 Markov Chain Monte Carlo search for best-fit models 



For each of the parametric models, there are at least six parameters to constrain. We use Bayesian inference to provide a 
probabilistic distribution of tri-axial bar model parameters by combining prior information with the data. Information from 
the data and prior knowledge entirely determines the posterior probability distribution as follows: 

n(x) oc £(data|a;) x prior(a;). (10) 

We use Markov Chain Monte Carlo (MCMC) simulations to explore the parameter space and find the posterior distri- 
bution. MCMC can be implemented in a nu mber of ways to explore the parameter space. In this work we use the standard 



random walk Metropolis-Hastings algorithm ( Metropolis et al. 19531 : Hastings! 1970T ). At any point x n in the chain we gener- 
ate a new trial point atrial by drawing parameter shifts from the symmetrical proposal distribution function <j>(x — x n ). This 
algorithm is universal for any proposal distribution function. Here, we use a multi-dimensional Gaussian distribution in the 
vector of parameters. We accept a candidate sta te x n +i with probability min{l, 7r(x tl -i a i)/7r(x n )}, where n(x) is the equilibrium 
distribution of the chain ( Lewis fc Bridle! 20021 ) . 



For each model, the chain is started at the corresponding best-fit model values of parameter vector as given in lRattenburv et al 



(20071). The model total star counts are normalised so that they reproduce the observed iVo for each sightline. The MCMC 
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Model 


Sightlines 


Ro 




2/0 


zo 


a 




SI 


8b 


xl 


xl 


xl 


*4 


X 2 /d.o.f 


Gi 


7631 


8.22 


1.23 


0.51 


0.43 


28.8 








13801 


15858 


15349 


2300 


2.07 


K 2 


8118 


8.17 


0.68 


0.28 


0.26 


31.7 








12384 


17710 


11714 


2691 


1.83 


E 3 


8271 


8.14 


0.69 


0.29 


0.28 


28.8 








12316 


18698 


11302 


937 


1.74 


G'l 


7749 


8.27 


1.24 


0.52 


0.42 


29.1 


-0.88 


-0.31 


-0.012 


13995 


16235 


15513 


2253 


2.07 


E' 2 


8178 


8.21 


0.68 


0.28 


0.25 


31.8 


-4.73 


-0.22 


-0.023 


12587 


17902 


11739 


2168 


1.81 


V 3 


8196 


8.20 


0.81 


0.33 


0.30 


32.6 


-7.12 


-0.21 


-0.018 


12825 


18031 


11314 


2228 


1.81 



Table 1. The best fit values of six- and nine-parameter tri-axial galactic bar models. 



simulation is then allowed to run, varying each of the six (or nine) parameters. Given the sampling strategy described above, 
we also use an adaptive step size to guarantee a sufficient sampling of the posterior, leading to convergence to the equilibrium 
posterior distribution in the chain after several thousand steps. We run 10 5 chains for each model to obtain the posterior 
distribution of the parameters, after neglecting the first half of the chain to allow the chain to reach its equilibrium state. 



4 RESULTS 

For the simple bar models, we have six parameters Ro,Nq,xo, yo, zo and a. We measure 2.94 x 10 6 RC stars in the OGLE-III 
data over a viewing area of 90.25 deg 2 in total. The total number of sightlines is 9019. For each model, we first run 10 5 Monte 
Carlo chains and get the best fit parameters for each model based on the original data set. We exclude sightlines that deviate 
more than 4cr, i.e. where (X™ od - X° bs ) T Sr 1 (X™ od - X obs ) ^ 16, from the model and then re-run to obtain the final results. 

The results of best fit parameters are presented in Table Q] Since No is not a true parameter and can be normalised by the 
total star counts, we include No in the MCMC chain but do not give its value here. The number of sightlines for each model 
is also listed in TabreQ] We separate \ 2 contributions from the number counts (xf), mean distance modulus (xl)j geomertical 
variance (xi) respectively. \\ is the contribution from non-diagonal terms in the error matrix. Among these models, the E2 
and E3 models have similar \ 2 /d.o.f values and they are both much smaller than that of the Gi model. This suggests that 
the Ei and E3 models best fit the data, while the Gi model fits the data worst. 

We present the marginalized posterior distribution of the parameters for the six-parameter Gi model in Fig. [1] We can see 
that the probability distribution function (PDF) of each parameter is nearly Gaussian. We also show the degeneracy between 
each parameter in Fig. [T] The 1- and 2-cr regions (i.e., 68% and 95% confidence levels) are shown in inner and outer contour 
lines respectively. Fig. [5] is the same as Fig. [T] , but for the E2 model. 

We compare the model prediction quantities, e.g. the surface density of stars T,rc, mean distance modulus fj., and 
geometrical variance cr 2 with that directly from observation for the models Gi and E2. In the top left panel of Fig. [2l 



we 



plot T,rc for three stripes at |6| = 2.5°, 4.0°, 5.5° as observed in Nataf et al. ( 20121 ). E_rc is maximized near I = 0° for all 



latitude stripes. In the same figure, we also plot the Gi model predictions of £_rc- At small longitudes, the model predictions 
are consistent with the observations. However, at large longitudes, especially for the |fe| = 5.5° stripe, the model predictions 
are systematically lower than what is observed. In the top right panel of Fig. [2l we plot T,rc from both model predictions 
and observations at fixed longitudes of / = 0.0°, 5.0° and 10.0° respectively. It shows the expected result that number counts 
increase monotonically with decreased separation from the plane for each longitude. However, it again shows that at large 
longitude, especially for the I = 10.0° stripe, the model predictions T,rc are systematically lower than the observed values. 

The middle left panel of Fig. [2] shows the Gi model (black symbols) and observed (red symbols) mean distance modulus 
fj, against longitude. In the middle right panel, we plot /i for both model prediction and observation at fixed longitudes of 
/ = —5.0°, 0° and 5.0°. In the bottom panels of Fig. [6] we plot the modelled (black symbols) and observed (red symbols) 
geometrical variance cr 2 vs. longitude (left) and latitude (right) respectiv ely. 

Notice that negative values of cr 2 , should not be a cause of concern. iNataf et alj (j2012T j estimated the dispersion due to 



Z Z Z Z /t t\ 

°" M = a -a I -a E (11) 

where a 2 is the measured brightness dispersion of the RC, 07 = 0.09 is the estimated intrinsic magnitude dispersion of RC, 
and a% is the dispersion due to differential extinction. Therefore, for sightlines where the measured brightness dispersion is 
small relative to the intrinsic dispersion and differential extinction, we can sometimes obtain negative variance in cr 2 . We do 
not remove these sightlines with negative geometrical variance to avoid systematically biasing the final output to have large 
variance. However, in the bottom panels of Figs. [2] and [6] the predicted cr 2 are well concentrated while the data have very 
wide variations. It could be that the models are wrong, or it could be that the data are wrong. To test for the latter, we 



take the intrinsic dispersion of the RC as a free parameter instead of 07 = 0.09 as quoted in INataf et al. (2012). Though very 
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Model E2 with cut 


Sightlines 


Ro 


.f<) 


yo 


ZQ 


a 


xl 


xl 


xl 


xl 


X 2 /d.o.f 


\l\ ^ 4° 
\l\ < 4°,|6| < 4° 


5439 
4033 


8.17 
8.18 


0.80 
0.71 


0.25 
0.24 


0.24 
0.26 


31.4 
30.8 


8127 
5432 


11525 
9085 


8166 
5472 


1039 
237 


1.77 
1.67 



Table 2. The best fit values of tri-axial galactic bar model after removing the outer part of the bulge regions. 



different values of crj are suggested for different models; the models show little improvement in the x 2 /d.o.f, and thus will 
not be discussed any further. 

Comparing Fig. |S]with Fig. [5] we can see that a clear improvement over the Gi model is that at large longitude, e.g. for 
the |6| = 5.5° stripe, the model surface density of stars E_rc is more co nsistent with the o bserved data. 

In Fig. [3] we plot the mean distance to RC centroids measured by Nataf et al. (2012) and give the G\ model predictions 
in the same figure for comparison. The distances are projected onto a face-on view of the central region of the Galaxy. The left 
panel of Fig. [3] shows data for sightlines with 2.5° ^ |b| ^ 3.5°, while the right panel is for sightlines with 4.5° ^ |6| ^ 5.5°. 
Fig. [7] is the same as Fig. [3] but for the E2 model. The E2 model is broadly consistent with the observations, but does not 
yield a very significant flattening at small positive longitudes as is seen in the data. 

An interesting phenomenon, re cently discovered, is that the RC of the bulge is observed to sho w two peaks in distance 



at large separations from the plane ( Nataf et al. 2010l ; McWilliam fc Zqccali 



2010; 



shown to be due to an X-shape structure by means of N-body models (|Li fc Shen 



Sait o et al 
2012| : 



2011 



Ness et al 



, a re sult that has been 
20121 1. The existence of 



the vertical X-shaped structure hints that the formation of the Milky Way bulge is shaped mainly by internal disk dynamical 
instabilities. 

To verify the X-shaped fine structure, we calculate the excess fraction in star counts between the model and the observed 
data, (Nobs — iV mo d e i)/Mnodei, for each sightline. We then project the excess fraction to produce a 2-D residual map, shown 
in the top panel of Fig. [4] for G\ model. We can see a dim X-shaped structure extending from I = 13° to I — —10°, with 
latitudes \b\ ^6°. After removing the bad sightlines, we still have 2.52 x 10 6 stars remaining. As we find a best-fit predicted 
RC population for the bulge of 2.35 x 10 6 , the excess RC population is 7.3%. In the bottom panel of Fig. 2] we also present 
the residual map of mean distance modulus. The regions have larger number counts show larger distances in observations 
than in models, but this trend does not hold on a pixel- by- pixel basis. Fig. [8] is the same as Fig. [4] but for the E2 model. 
The predicted RC population for the bulge is 2 .35 x 10 6 while the observed value is 2.52 x 10 6 with an excess of 5.3%. This 
value is not inconsistent with Li fc Shen |2012l ). who found the excess light to be 7% in simulation. The agreement may be 
somewhat coincidental since their values are for an edge-on view. The existence of an X-shape structure suggests that at the 
outer regions of the bulge, there is some fine structure in the Galactic bar. 

We also test the results of these models after introducing a centroid offset (51, 5b) and tilt angle (/3). To avoid confusion, 
these models are labelled as Gi,E 2 and E' 3 . The results of the nine-parameter model fitting are also shown in Table [1] 
However, for the G'\ and E' 2 model, the results of fitting a nine-parameter model show little improvement over the simpler 
six-parameter models. For the nine-parameter E' 3 model, a very slightly larger \ 2 /d.o.f is obtained due to the more sightlines 
being included. This is perhaps not surprising since our sightlines do not cover the very central regions which presumably 
show the most prominent deviations from a non-offset model. Thus from now on /3, 51 and 5b will all be fixed at zero. 



4.1 Disk contamination 



Ther e is evidence that there is some structural change in the Galactic bulge inward of I ^ 4° (jGerhard fc Martinez- Valpuesta 



20121 ). As mentioned above, the vertical X-shaped structure is assumed to be shaped mainly by internal disk dynamical 



instabilities and exists at larger longitudes. Both features are a source of disk contamination which is difficult for our model to 
reproduce, especially in regions remote from the centre of the bulge. From Figs. [2] and [6] we can see that at larger longitude, 
the number counts and geometrical variance in observation also systematically deviate from the model prediction. 

To see the effects of possible disk contamination, we simply exclude data from the outer part of the bulge. We adopt two 
cut criteria. Firstly, we remove sightlines with \l\ > 4°. Secondly, we retain only the sightlines whose longitude and latitude 
are both smaller than 4° . The value of 4 degrees is chosen because in this region, there is a slope change in the particle surface 
density, see, e.g., Fig. 2 of iGerhard fc Martinez- Valpuesta I (|2012h . 



We test the results of the E2 model for the above two cut criteria. The results are given in Table [2] The model fit has 
an obviously smaller \ 2 /d.o.f than the fit without cut. We again calculate the model-predicted RC population, and find the 
excess reduces to 3.3% for cutting the area of |2| > 4° and 2.6% for cutting |Z| > 4° and |6| > 4°. These are not surprise since 
we can see from Figs. [6] and [8] that the larger dispersion and excess both lie in the outlier of the bugle region. It suggests that 
we can get a much more realistic inner bar structure after cutting the wide longitude sightlines. The cuts on the data partly 
avoid the disk contamination and fine structure, which are presently very difficult to handle. 
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4.2 Comparison with results from OGLE-II 



Rattenbury et al.l (|2007i ) used RCG stars from 44 bulge fields from the OGLE-II microlensing collaboration database to 
constrain analytical tri-axial models for the Galactic bar. They found the bar major axis is oriented at an angle of 24° to 27° 
to the Sun-Galactic centre line of sight. They also found the ratio of semi-major and semi-minor bar axis scale lengths in the 

Galactic plane xo, yp, and vertical bar scale length zp, to be xp : yp : zp = 1.00 : 0.35 : 0.26. 

As mentioned in the introduction, the data from OGLE-III as analyzed by Nataf et al. ( 20121 ') offer many more sightlines 
and better calibration in the RCG population. Both are helpful to refine the constraints on the bar parameters. We remove 
the area in the outer part of the bulge region and keep the fields with —4° ^ I ^ 4° and —4° ^ b ^ 4° to partly remove the 
disk contamination and concentrate on the inner region of the bar. When using all the data, we found an excess, consistent 
with the so-calle d the X-shape structure in the outer part of the bulge region. The excess is ~ 7%, in agreement with 
Li fc Shenl |2012h , but which complicated the fitting. When we remove data from the outer part of the bulge, the effect of 



the fine structure is reduced and we can fit a smoother inner bar structure. We obtain a bar angle oriented at 29° to 32°, 
and the axis ratios are around xp : y p : zp = 1.00 : 0.41 : 0.38. This suggests a bulge that is not as tri- axial as the bar 
structure found by Rattenbury et al. (2007); it is also vertically more extended than the working model of Gerhard ( 20021 ) 
with xp : yp : zp = 1.00 : 0.40 : 0.30. 



5 DISCUSSION 



We use the RCG population in the OGLE-III database to fit a tri-axial Galactic bar model and to derive much tighter 
constraints on the bar parameters. When using all the available data, we found there is evidence fo r an X -shaped structure in 
the outer part of the bulge region, with observable properties consistent with those of Li fc Shen ( 20121 ). This fine structure 



complicated the fitting of our analytic models. We removed data from the area in the outer part of the bulge region and 
retained the fields with —4° ^ I ^ 4° and —4° ^ b ^ 4° to partly remove disk contamination and to concentrate on the inner 
region of the bar. This reduced the effect of the fine structure and allowed a fit to the resulting smoother inner bar structure. 
We obtain a bar angle between 29° to 32° with an axis ratio a round xp : yp : zp — 1.00 : 0.41 : 0.38. This bar is close to 



being prolate, and no t as triax i al as the bar structure found by iRattenburv et al.l (|2007l ); it appears vertically thicker than 
the working mode l of Gerhard ( 20021 ) who gives xp : yp : zp — 1.00 : 0.40 : 0.30. The bar is more prolate than that found by 



Dwek et al. 1 1 19951 ) and has a larger bar angle. The m odel derived here can effectively reproduce the number density, distance 
modulus over the sky probed by the OGLE-III data 1 Nataf et al. 20121 ). 

Future microlensing observations will provide us with additional constraints on the structure of the Galactic bar. Based 
on the improved and tightly constrained bar model, it will be interesting to re-calculate the optical depths and time-scale 
distribution for microlensing events, and see whether they are consistent with the data. 

The X-shaped structure appears in both simulations and in our data. Its existence in the Milky way implies that the 
Galactic bulge is shaped mainly by internal disk dynamical instabilities instead of mergers (e.g. iLi fc Shenl |2012| ). How to 
model subtle structures such as the X-shaped excess remains a challenging task. One possible way forward is to expand the 
density di stribution into ortho gonal basis functions; a similar approach is taken in weak lensing to model complex point spread 



functions (jAndrae et al.ll201ll ): we plan to explore this in the near term. 
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Figure 1. The marginalised probability density distributions (PDFs) and contours in the plane of Ro, xo, yo,zo, a for the six-parameter 
Gi model. All the lengths are in units of kpc and the angle a is measured in degrees. The inner and outer contours are respectively for 
the 1-tr and 2-cr confidence level regions (for two parameters simultaneously). 
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Figure 2. The predicted surface density of stars mean modulus [i, and geometrical variance cr^, vs. those directly from observation 

Gi model. Top panels: the surface number density T,rq vs - longitude (left) and latitude (right). Black symbols denote model values 
while colour symbols denote observational data. Middle panels: mean distance modulus fi vs. longitude (left) and latitude (right). In the 
middle and bottom panels, black symbols denote model values and red symbols observational data. Bottom panels: geometrical variance 
a 1 ^ vs. longitude (left) and latitude (right). 
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Figure 3. The face-on view projection of the distances to RC centroids for the Gl model. The red points show the projected location 
of each measured RC centroid, and black points show the model prediction. The left panel is for sightlines satisfying 2.5° ^ |6| ^ 3.5°, 
while the right panel is for 4.5° ^ |b| ^ 5.5° sightlines. 
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Figure 4. The residual map of number counts (top) and mean distance modulus (bottom) for the six-parameter G\ model. The residual 
of number counts is defined as the fractional difference in star counts between the model and the observed data, (-/V t> s — JV mo( i c i) /iV mo del- 
The total excess fraction of number counts is 7.3%. The residual in the mean distance is defined as £i b B — Mmodel = 5 log ^obs/^modeb 
where D b s and -D mo dcl are the the observed and predicted distances. 
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Figure 6. The same as Fig. [2] but for the six-parameter E2 model. 
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Figure 7. The same as Fig. [3] but for the six-parameter E2 model. 
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Figure 8. The same as Fig. [4] but for the six-parameter E2 model. The total excess fraction of number counts is 5.3%. 
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